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COPYRIGHT («) 1978, 1980, 1982, 1984 BY « 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. * 
ALL RIGHTS RESERVED. * 
® 

THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED * 
ON IN ACCORDANCE E AND WITH THE * 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER * 
AY NOT BE PROV ERW I . 

OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
TRANSFERRED. ' 
® 
THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
CORPORATION. * 
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; FACILITY: BASIC Language Support 


; ABSTRACT: 


This module accepts two scaled packed decimal values and returns the 
correctly rounded (when required) result. This is done by co computing 
a 62 digit result and than extracting the appropriate digits 


: ENVIRONMENT: AST Reentrant 
; AUTHOR: Bob Hanek, CREATION DATE: 19-JAN-1981 
; MODIFIED BY: 


: - 01 - Original. RNH 01-JAN-1981 

; 1-002 - Change name to BASSEXTEND MULP and change ‘errors to 
Decimal Error. PLL 12-Feb-1982 

1-003 - Add code to check Ry 37 flags field within BASIC frame 
for decimal overflow detection, and feed that info into 
the PSW. LB 14-May-1982 
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pac 6-SEP-1984 BASRTL.SRC JBASEMULP.MAR; 1 
-SBTTL BASSEXTEND_MULP = Extended precision packed multiply 


+e 
FUNCTIONAL DESCRIPTION: 


This module accepts two scaled packed decimal values and returns the 
correc rounded (when required) result. si 

a 62 digit result and than extracting the appropriate digit S. is 
particular, given two_ packed decimal yorves - end B' of know tength 


s done by comput tng 


we extend them to. two ae t volves of n let A = A1*10*1 
+ AQ and 4 T1098" + Be oyero*ts + BO, yhere 3 
=< AO, 60, B10 os 10°15 - 1, hal Al =< 10°16 = 1. and O=< B11 =< 9. 
Then 
AtB = (A1#10°15 + AO) *(B1#10°15 + BO) 

= A1*B1*10*30 + (A1*B0 + hy Hata SB + AQ*B0 

= A1l*(B11#10*15 + B810)*#10*30 + (A1*BO + AOQ*B1)*10°15 + AO*BO 

= A1*B11*#10*°45 + A1*B10*10*°30 + (A1*BO + AQ*B1)*10°15 + AO*B0. 
Note that all_of the products of the form An*Bm have absolute value 
ha ey 10°31, so that each product can be exactly computed by a MULP 
nstruction. 


CALLING SEQUENCE: 


CALL BASSEXTEND_MULP(A.rp.dsd, B.rp.dsd, C.mp.dsd ,RND_TRUNC.rb.v) 


FORMAL PARAMETERS: 


A and B are the scaled decimal source strings 

C is the scaled decimal result string 

RND_TRUNC is a flag Indices ine, whether the final answer should be 
Founded or truncated. RND_TRUNC = 5 gives rounded result and 
RND_TRUNC = 0 gives truncated result. 


IMPLICIT INPUTS: 


NONE 


IMPLICIT OUTPUTS: 


NONE 


SIDE EFFECTS: 


Signals PACKED OVERFLOW whenever the product is too large to fit in 
the scaled decimal result cperand 
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08 BASSEXTEND PULP - extended precision pac 6-SEP-1984 99:38:39 BASRTL.SRCJBASEMULP.MAR; 1 ee 
O3FC 145 -ENTRY BASSEXTEND_MULP, “M<R2, R3, RG, RS, R6, R7, RB, RI> 
| 139 3 Entry point 
| re ; Move parameters to registers R6 through R9 
2 AC 7D 09 139 MOva 4(AP), R6 ; Move first two parameters to R6/R7 
| Cac 7D 44 123 MOVa 12(APS, RB ; Move next two parameter to R8/R9 
BRA 138 ; REGISTER USAGE: 
| OOA 155; RO Through R5 are reserved for the packed decimal instruction 
OOA 136 : R6: Pointer to A descripter 
OOA 157 ; R?7: Pointer to B descripter 
OOA 158; R8: Pointer to C descripter 
OOA 159; R9: Round/trun flag 
44 9 160 ; 
O0OA 161 
QO0A 166 H 
QOOA 1635 ; Allocate stack storage and turn off decimal overflow reporting. 
OOA 164 ; Check the setting of the decimal overflow bit in the flags field 
OOA 165 ; within the BASIC frame. Use that value to set the PSL decimal 
4 9 166 ; overflow trap setting (enable or disable). 
5E 0000007C 8F ce QOOOA 168 SUBL #TOT, SP ; Allocate 26 longwords on stack 
54 OC AD DO 0011 169 MOVL 12(FP) RG ; Fetch the saved frame pointer 
55 00000000°'GF DE 0015 170 MOVAL G*“BASSHANDLER,RS ; Fetch addr of BASIC handler 
55 64 D1 OO1C 171 CMPL O(CR4) RS ; Check if this is a BASIC frame 
06 13 OO1F 78 BEQL 5$ ; Branch if a BASIC frame 
0080 8F 68 0021 17 BISPSW #*x80 i; Else, set DV bit in PSL 
17 11 0025 174 BRB 15$ 
54 £6 A4 B80 0027 175 S$: MOVW =26(R4) RS : Fetch the flags word 
54 FBFF of AA 0028 178 BICW #*XFBFF RG ; Mask out all but DV bit 
4 B65 00350 17 TSTwW RG ; Check if DV was set in BASIC frame 
06 13 0032 178 BEQL 10$ ; Br if not enabled 
0080 8F B88 0034 179 BISPSW #*x80 ; Set DV bit in PSL (logical OR) 
11 0038 180 BRB 15$ ; Continue as usual 
0080 8F B89 OO3SA 181 10$: BICPSW #*x80 ; Clear DV bit in PSL (logical AND) 
78 AE OC O03E 13 15$: MOVPSL DVFLAG(SP) ; Save current PSL 
78 AE FF7F 8F AA 0041 18 BICW @*XFF7F, DVFLAG(SP) : Save only decimal overflow bit 
0080 8F B9 Shee is BICPSW #*x80 ; Turn off decimal overflow reporting 
Ohee 136 ; Break B into its high and low parts 
6E 1F 00 04 867 67 OO F8 pe8 138 ASHP #0, (R7), A4(R7), #0, #31, (SP) ; ertons B to 31 gigits 
70 AE «OB AE «67D «(0053—18 mova 8 §P), gb«sP) ; Set BO to low 13 digits of B 
10 OO 6€ 6 ¢€F18F £8 0058 «190 ASHP #-15, #31, (SP), #0, #16, BI(SP); Set Bl to high 16 digits of B 
64 AE OO5F 
061 191; , 
bet 136 ; Break A into its high and low parts 
6€ 1F 00 04 86 66 00 F8 0061 194 ° ASHP #0, (R6), O4(R6), #0, #31, (SP) ; Extend A to 31 gigits 
SC AE OB AE 7D 0069 195 mova 8($P), Ad(SP) ; Set AO to low 15 digits of A 
10 00 6E€ IF 4 4 F8 oF 196 ASHP #-15, #31, (SP), #0, #16, AI(SP); Set Al to high 16 digits of A 
7 197 ; 
305 ; Compute A0*BO and store in the low digits of the result. 
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1 7O AE OF SC AE +4 25 200 MULP #15, AOCSP), #15, BO:SP), #31, LOW(SP) 


; Compute A0*B1 and store the results in a 31 digit temporary 
: MULP =s-« #15, AO(SP), #16, BI(SP), #31, TI(SP) 


mrorore 


Fas 


1f 64 AE 10 SC AE OF 25 


mmtetatentesn=nen in 


20 AE 
0 > 3 
| 0 Og g ‘ d low 16 Gigits of A0*B1 to low digits of the result and move the high 
0 6 i 1 digits of A0*B1 to high digits of the result 
| fF 00 20 AE af +4 F8 4 § . ASHP #15, #31, TICSP), #0, #31, T2(SP) ; T2 = Lo 16 digits of A0*B1 
| 1F 00 20 AE 1F Hi fl F8 210 ASHP #-16, #31, TICSP), #0, #31, (SP) ; (SP) = Hi 15 digits of A0*B1 
10 AE 1F 30 AE fF 20 0099 11 ADDP #31, T2(SP), #31, LOW(SP) ; Add lo 16 digits to lo 
06 1C QOA4 \¢ BvC 1$ ; digits of the result and 
on ce § FO ee Bae F ACOP #1, MONE, #31, (SP) i and propagate the carry 
DoAt 15 : Compute A1*B0 
1F 70 AE OF SO AE 10 25 OAC 3 is: MULP =—s«#'16, A1(SP), #15, BO(SP), #31, TI(SP) 
de i§ : . d low 16 qigits of Al*B0 to low digits of the result and add the high 
s+ 3! 3 digits of Al*B0 to high digits of the result 
1F 00 20 AE 1F 9 F8 ie44 22 ’ ASHP #15, #31, TICSP), #0, #31, T2(SP) ; T2 = lo 16 digits of Al*B0 
iF 00 20 AE fF f° 4 F8 a4 223 ASHP #-16, #31, TI(SP), #0, #31, T3(SP); T3 = hi 15 digits of Al*B0 
10 AE IF SO AE 1F 20 QOC9 224 ADOP #31, T2(SP), #31, LOW(SP) ; Add lo digits to lo digits 
07 iC 00D0 5 BVC : Branch if no carry 
40 AE 1F 1C 8F 01 $ Ode $ ADDP #1, #ONE, #31, T3(SP) ; Propagate carry 
F 4OAE 1F 20 9008 ? 2$: ADDP Oss #34, T3(SP), #31, (SP) ; Add hi digits to hi digits 
oF > Compute 811*10*15 and A1*B10 
40 AE 65 AE OF 34 OO0DF 7 movP a5. B10(SP), T3(SP) ; 13 = 810 
64 AE 10 40 AE OF $ 0E5 3 SUBP #16, BI(SP) Bl = 811*10°15 
1F 50 AE 10 40 AE OF OEC MULP a8" I (SP). #16, AI(SP), #31, 118) ; T1 = A1*B10 
Sal Bese 
F6 35 3; Add tou ¢ digit of A1*B10 to low digits of result and add high 30 digits of 
ae $ 3; Al*B1 igh digits of result. 
| 1F 00 20 A€ of 1E F8 OF6 3 : ASHP =: #30, «#31, TI(SP), #0, #31, T2(SP) ; T2 = lo digit of A1*B10 
fF 00 20 AE 1F AA S F8 o107 239 ASHP #-1, #31, TICSP), #0, #31, T3(SP) ; T3 = hi 30 digits of Al*B10 
10 AE 1 30 AE fF 40 109 40 ADOP #31, T2(SP), #31, LOW(SP) ; Add lo digit to lo digits 
$6 C 0110 41 BvC 3 ; Branch if no carry 
6 fF 1C 8F 1 $0 \1¢ 4 ADDP al #31, (SP) : aropagete ¢a Ay aH 
6€ IF 40 1F 20 118 43 3$: ADDP #31, itty #31, (SP) : Add hi digits to hi digits 
NE $2 : Compute A1*B11*10°14 and add to high digits of result 
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OF 00 646 AE 10 44 +4 F8 1) 247 ASHP #-1, #16, BI(SP), #0, #15, T3(SP) ; T3 = 811#10*14 
| 50 AE 10 40 Ato +4 25 } 248 MULP #15, T3(SP), #16, AI(SP), #31, TICSP) ; T1 = A1#B11¢10°14 
| 
| 6E fF 20 AE 1F 20 : rH ADDF #31, TICSP), #31, (SP) ; Add to high digits 
| : af : Determine the number of place we need to shift the result. 
55 08 A? 4 “6 8) 1 34 . AdDB3 pcre), 8(R7), RS ; RS = total # of fraction digits 
08 A 2 iF 22 SUBB (RB), RS ; RS = # of extra fractional digits 
| 14 2§ ; Shift the result the necessary number of digits, round if necessary, and 
12 2 ; Store the result in the destination operand. 
fF 59 10 AE ¥ oo 3 F8 Bie 39 : ASHP RS, #31, LOW(SP), RO, #31, TI(SP) ; T1 = final result lo digits 
55 1F 80 148 60 ADDB #31, R5 3; RS = shift for high otetts 
30 AE 1F 00 =«=6E iF 55 F8 14E 61 ASHP ae #31, (SP), #0, #31, T2(SP) ; T2 = final result hi digits 
30 AE «O1F «6020 AE OTF «020s 0156 6¢ ADpP = ss w#34, TI(sP), #31, °T2(SP) : 12 = Final result 
16 1D 015D BVS OVERFLOW : Branch if overflow 
68 00 30 AE IF + F8 HEA 64 ASHP #0, #51, T2(SP), #0, (RB), a4(RB) ; Store result in destination 
0B 1D 168 65 BVS OVERFLOW ; Branch if overflow 
AE B8 OQ16A ¢ RETURN: BISPSW DVFLAG(SP) ; Restore decimal overflow flag 
5E 0000007C BF CO 190 6 ADDL #TOT, SP ; Deallocate stack storage 
04 0174 68 RET ; Return 
B172 69 ; 
017 70 ; 
0175 71; 
0175 272 OVERFLOW: 
78 AE 685) «(0175 7 TSTW DVFLAG(SP) ; Did caller enable Dec ovfl reporting? 
FO Fs 178 74 BEQL RETURN :; Branch if reporting not enabled 
78 AE 8B 17A 75 BISPSW DVFLAG(SP) ; Restore decimal overflow flag 
5E 7¢ or CO 017D 7 ADDL #TOoT, SP ; Deallocate stack storage 
00000000'8F DD 0184 27 PUSHL #BAS$K_DECERR : Signal overflow 
100 * GF 01 FB 018A 78 CALLS #1, G*BASS$S$STOP 3 
04 0191 279 RET : 
019 80 
019 81 
019 82 END :; End of module BASSEXTEND_MULP 
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B10 

BASSSSTOP eeeerere x 0 

BASSEXTEND MULP 00000000 RG 1 

BASSHANDLER etereree x 0 
00 


@ 
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BASSK DECERR eeereree x 
DVF LAG 

00001 
ONE 

OVERFLOW 175 R 4 
_ O0¢0 R 1 
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18 
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eee cen ewes eases + 


PSECT name Allocation PSECT No. Attributes 


. ABS 2 
-BASSCODE 00000192 (« 402.) O01 ¢ 1.) PIC USR CON REL LCL 


Phase Page faults CPU Time Elapsed Time 
Initialization 0 00:00:00.07 0:00:00.30 
Command processing 120 0:00:00.50 0:00:01.83 
Pass 5 BS BN O88 ae 
Symbol table sort ! Bp: 7:00.01 0:00:00.01 
Pass 6 BOR: 3-83 00:00:01.17 
Symbol table output 00:00:00.02 Bo + 80 2 80-08 
Psect synopsis output 00:00:00.01 00:00:00.01 
Cross-reference output 00:00:00.00 00:00:00.00 
Assembler run totals 299 00:00:02.16 00:00:05.99 


The working set Limit was 750 pages. 

5342 bytes (11 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 18 non-local and 6 local symbols. 
282 source Lines were read in Pass 1, producing 11 object records in Pass 2. 

0 pages of virtual memory were used to define 0 macros. 


AX/VMS Macro V04-00 
BASRTL.SRCIBASEMULP AR; 1 


00000000 < 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
E RO NOWRT NOVEC LONG 
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varelt Macro Run Statistics 984 :29 [BASRTL.SRCIBASEMULP.MAR; 1 (4) 


temo coooeeocoeooseoesococes$ 


! Macro Library statistics ! 


tent r cee seme moe were rere wero ee + 


Macro Library name Macros defined 


-$255$DUA28: CSYSLIBJSTARLET.MLB; 2 0 
O GETS were required to define 0 macros. 


There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESS1ON/D1 SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:BASEMULP/0BJ=0BJ$:BASEMULP MSRC$:BASEMULP/UPDATE=(ENH$:BASEMULP) 
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